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This paper describes the performance and sensitivity to neutrino mixing parame- 
ters of a Magnetised Iron Neutrino Detector (MIND) at a Neutrino Factory with a 
neutrino beam created from the decay of 10 GeV muons. Specifically, it is concerned 
with the ability of such a detector to detect muons of the opposite sign to those 
stored (wrong-sign muons) while suppressing contamination of the signal from the 
interactions of other neutrino species in the beam. A new more realistic simulation 
and analysis, which improves the efficiency of this detector at low energies, has been 
developed using the GENIE neutrino event generator and the GEANT4 simulation 
toolkit. Low energy neutrino events down to 1 GeV were selected, while reducing 
backgrounds to the 10~ 4 level. Signal efficiency plateaus of ~60% for and ~70% 
for events were achieved starting at ~5 GeV. Contamination from the — > u T 
oscillation channel was studied for the first time and was found to be at the level 
between 1% and 4%. Full response matrices are supplied for all the signal and back- 
ground channels from 1 GeV to 10 GeV. The sensitivity of an experiment involving 
a MIND detector of 100 ktonnes at 2000 km from the Neutrino Factory is calculated 
for the case of sin 2 20 is ~ 10 _1 . For this value of #13, the accuracy in the measure- 
ment of the CP violating phase is estimated to be A5cp ~ 3° — 5°, depending on 
the value of 5c p, the CP coverage at 5<7 is 85% and the mass hierarchy would be 
determined with better than 5a level for all values of Sep- 
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I. INTRODUCTION 

The Neutrino Factory, a new type of accelerator facility in which a neutrino beam is 
created from the decay of muons in flight in a storage ring, is perhaps the most promising 
facility design to resolve the problem of CP violation in the neutrino sector. The physics 
potential of this facility was first described by Geer [T]. The expected absolute flux and 
spectrum of neutrinos from such a facility can be calculated with smaller systematic errors 
than those associated with the beams of alternate facilities due to the ability to measure the 
muon beam flux and the highly accurate measurement of muon decay kinematics [2]. Since, 
in principle, both and \i~ can be created with the same systematic uncertainties on the 
flux, any oscillation channel can be studied with both neutrinos and antineutrinos, improving 
sensitivity to CP violation. Table [I] shows the oscillation channels that will contribute to 
the flux at any far site due to the decay of 

TABLE I: Oscillation channels contributing to flux from the decay of fi + . 



v e origin 



origin 



v e — > v e (u e disappearance channel) 
v e v fi (Golden channel) 



v r (Silver channel) 



(v^ disappearance channel) 
— > V T (Dominant oscillation) 



—7- v e (Platinum channel) 



The sub-dominant v e — > oscillation [3] was identified as the most promising channel 
to explore CP violation at a Neutrino Factory. The charged current interactions of the 
"Golden Channel" produce muons of the opposite charge to those stored in the storage 
ring (wrong-sign muons) and these can be detected with a large magnetised iron detector [3] . 
The original analyses were carried out assuming a Neutrino Factory storing 50 GeV muons 
and, as such, were optimised for high energy using a detector with 4 cm thick iron plates 
and 1 cm scintillator planes. However, subsequent phenomenological studies carried out as 
part of the International Scoping Study (ISS) for future neutrino facilities [HI E] favoured 
a stored muon energy of 25 GeV and showed the importance of neutrinos with energies 
below 5 GeV. The Magnetised Iron Neutrino Detector (MIND) is a large scale iron and 
scintillator sampling calorimeter, similar to MINOS [7], which was re-optimized from the 
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original studies motivated by these findings [8HT0]. The performance obtained indicated 
that the combination of two Magnetised Iron Neutrino Detectors at 4000 km and 7500 km 
would give optimum sensitivity to the mixing parameters [TT] . 

The studies of MIND mentioned above evaluated the performance of the detector us- 
ing deep inelastic scattering events only, with a simplified simulation, reconstruction and 
kinematic analysis. The performance needed to be evaluated and improved using a full sim- 
ulation and analysis of all physical processes. As part of the International Design Study for 
a Neutrino Factory [121 03] , a software framework to perform these studies has been devel- 
oped. Pattern recognition and analysis algorithms were developed and first applied to data 
generated using the same simulation as was used in the ISS studies. The development of 
the algorithm and the results of its application were described in [T4] and [15], where it was 
shown that under these conditions the efficiency and background could be maintained at a 
similar level to that achieved in the ISS studies. This paper introduces the full spectrum of 
possible neutrino interactions generated using the neutrino event generator GENIE [16J and 
a comparison with another event generator, NUANCE [T7] . These interactions were tracked 
through a new GEANT4 simulation [181 HH] with full hadron shower development and a new 
detector digitisation not present in previous studies. The events were then subject to the 
pattern recognition algorithm presented in [15] . reoptimised for the new simulation. Finally 
a likelihood based analysis was used to further suppress backgrounds. A preliminary ver- 
sion of this analysis using the NUANCE package has been published in the Interim Design 
Report of the IDS-NF [20] . This paper includes the full GENIE simulation, a comparison to 
NUANCE, an estimate of systematic errors, and sensitivity calculations for #13, the neutrino 
mass hierarchy (sign of Am 2 3 = ml — tttJ) and the CP violating phase Sep- 

Recent results from the reactor experiments Daya Bay, RENO and Double Chooz [2TH23] , 
as well as evidence from T2K [21] and MINOS [25], have demonstrated that the value of 
#13 is large (with a combined average of sin 2 2^13 = 0.097 ± 0.012). These results increase 
the likelihood of a discovery of CP violation and the determination of the mass hierarchy in 
neutrinos. It was shown in the Interim Design Report of the IDS-NF J2D] that at a value of 
sin 2 2#i3 ~ 0.1 the optimum Neutrino Factory configuration is achieved with a muon energy 
of 10 GeV and with a far detector at a distance of 2000 km. While the Neutrino Factory 
was designed to discover CP violation for a large range of values of #13 (down to values of 
sin 2 26* 13 ~ 10 -4 ), it will be shown in this paper that it also offers the best chance to discover 
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CP violation and the mass hierarchy at large values of #13, regardless of whether Am^ is 
positive or negative (inverted or normal mass hierarchy). 

This paper is organised as follows. Section [IT] introduces the relevant backgrounds and 



contaminations of the golden signal and describes the required suppressions. Section III 



describes the simulation tools and gives a description of MIND and the assumptions still 



made for this study. The analysis is described in Section IV, with a detailed demonstration 
of all variables and functions used to identify signal from background. The results from this 
analysis, including signal efficiencies, background rejection capabilities and performance to 
the Up — > v T oscillation signal, are presented in Section [Vl A discussion of some of the 
systematic errors of the analysis is described in Section VI Finally, full sensitivity to #13, 



5c p and the neutrino mass hierarchy will be presented in Section VIII Response (migration) 
matrices of this detector system for all signal and background will be shown in the Appendix. 



II. SOURCES OF IMPURITY IN THE GOLDEN SAMPLE 

The primary sources of background to the wrong-sign muon search come from the Charged 
Current (CC) and Neutral Current (NC) interactions of the non-oscillating neutrinos present 
in the beam. Specifically, the CC interactions of (z7 M ) being reconstructed as (z/ M ), NC 
from all neutrino types in the beam being reconsructed as [v^] and the CC interactions of 
v e (y e ) being reconstructed as (u^). Since these interactions are in far greater abundance 
than those of the signal channel and contain little or no discernible information about the 
key parameters # 13 and 5c p, they must be suppressed sufficiently so that the statistical 
error on the background is smaller than the expected signal level. This corresponds to a 
suppression of at least 10~ 3 for each channel in the signal region. 

In addition to the golden channel appearance oscillation there are three other appearance 
channels that will introduce neutrinos to the flux incident on the far detectors (shown 
in table [T]). The dominant oscillation, v^ip^) — > v T (v T ), which must be considered when 
fitting for the (v^) disappearance signal, should not pose a problem for fitting the golden 
channel since, at large #13, the dependence of this channel on this mixing angle is very 
small and the interaction would have to be reconstructed with the opposite charge to that 
of the true primary lepton. See (26] for a detailed discussion of tau contamination in the 
disappearance channel. The platinum channel, v^ip^) — > u e (P e ), should pose no problem 
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since the number of interactions should be similar to that produced by the golden channel 
and v e interactions produce a penetrating muon-like track in only a small fraction of cases. 
However, the silver channel oscillation, ^ e (^ e ) — > v T {v T \ would be expected to contribute a 
similar amount of v T to the flux as the golden channel does to u^, and since the primary r 
decays with a ~ 17.65% probablility via channels containing muons, a significant proportion 
of these interactions would be expected to pass the analysis cuts. As discussed in [27], 
fitting the observed spectrum without accounting for the presence of these v T interactions 
leads to significantly reduced accuracy in the fits. However, since this oscillation contains 
complimentary information about both 13 and Sep, handling this 'contamination' correctly 
has the potential to perhaps improve the fit accuracy compared to using an analysis which 
attempts to remove it from the sample. 

III. SIMULATION AND RECONSTRUCTION OF MIND 

In previous studies [Il[l0j[T5], only deep inelastic scattering (DIS) events generated with 
LEPTO 6.1 [28] were considered. However, at energies below 5 GeV there are large contri- 
butions from quasi-elastic (QE), single pion production (l7r) and other resonant production 
(RES) events. QE and In events are expected to exhibit lower multiplicity in the detector 
output, which makes candidate muon reconstruction simpler. This should improve recon- 
struction efficiency in low energy CC interactions but also potentially increase low energy 
backgrounds, particularly from NC In interactions. Other nuclear resonant events, produc- 
ing two or three pions, as well as diffractive and coherent production, have much smaller 
contributions. Moreover, the presence of QE interactions allows for the calculation of neu- 
trino energy without hadron shower reconstruction, improving neutrino energy resolution. 

A. Neutrino event generation and detector simulation 

Generation of all types of interaction was performed using the GENIE framework [T7] . 
The exclusive event samples generated by GENIE are shown in figure [TJ where 'other' 
interactions include the resonant, coherent and diffractive processes other than single pion 
production. The relative rates below 1 GeV are included for completeness, but will have 
negligible effect at a Neutrino Factory. GENIE also includes a treatment to simulate the 
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effect of re-interaction within the participant nucleon, which is particularly important for 
low energy interactions in high-Z targets such as iron. 

A new simulation of MIND using the GEANT4 toolkit [19J (G4MIND) was developed 
to provide flexibility to the definition of the geometry, to carry out full hadron shower 
development and to perform a proper digitisation of the events. This allows optimisation 
of all aspects of the detector, such as the dimensions and spacing of all scintillator and iron 
pieces, external dimensions of the detector and detector readout considerations. 

The detector transverse dimensions (x and y axes) and length in the beam direction [z 
axis), transverse to the detector face, are controlled from a parameter file. A fiducial cross 
section of 14 mxl4 m, including 3 cm of iron for every 2 cm of polystyrene extruded plastic 
scintillator (1 cm of scintillator per view), was assumed. A constant magnetic field of 1 T 
is oriented in the positive y direction throughout the detector volume. Events generated for 
iron and scintillator nuclei are selected according to their relative weights in the detector and 
the resultant particles are tracked from a vertex randomly positioned in three dimensions 
within a randomly selected piece of the appropriate material. Physics processes are modelled 
using the QGSPJ3ERT physics lists provided by GEANT4 [29]. 

Secondary particles are required to travel at least 30 mm from their production point 
or to cross a material boundary between the detector sub-volumes to have their trajectory 
fully tracked. Generally, particles are only tracked down to a kinetic energy of 100 MeV. 
However, gammas and muons are excluded from this cut. The end-point of a muon track is 
important for muon pattern recognition. 

A simplified digitisation model was considered for this simulation. Two-dimensional 
boxes - termed voxels - represent view-matched x and y readout positions. Any deposit 
which falls within a voxel has its energy deposit added to the voxel total raw energy deposit. 
The thickness of two centimetres of scintillator per plane assumes 1 cm per view. Voxels 
with edge lengths of 3.5 cm were chosen to match the required point resolution of 1 cm 
(3.5/\/l2), assuming a uniform hit distribution along the width of the scintillator bar. The 
response of the scintillator bars is derived from the raw energy deposit in each voxel, read 
out using wavelength shifting (WLS) fibres with attenuation length A = 5 m, as reported 
by the MINERvA collaboration [30]. Assuming that approximately half of the energy will 
come from each view, the deposit is halved and the remaining energy at each edge in x and 
y is calculated. This energy is then smeared according to a Gaussian with a/E = 6% to 
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FIG. 1: Proportion of total number of interactions of different v interaction processes for 
events generated using GENIE and passed to the G4MIND simulation. 
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FIG. 2: The digitisation and voxel clustering of an example event: (top left) bending plane 
view, (top right) non bending plane, (bottom) an individual scintillator plane. The 
individual hits are small dots (in red), the blue squares are the voxels and the black 
asterisks represent the centroid positions of the clusters. 



represent the response of the electronics and then recombined into x, y and total = x + y 
energy deposit per voxel. An output wavelength of 525 nm, a photo-detector quantum 
efficiency of ~30% and a threshold of 4.7 photo electrons (pe) per view (as in MINOS [7j) 
were assumed. Any voxel in which the two views do not make this threshold is cut. If only 



one view is above threshold, then only the view below the cut is excluded (see section III B ). 
The digitisation of an example event is shown in figure [2j 



B. Event reconstruction 

The reconstruction package was described in detail in j'fo i. We present here an update of 
the reconstruction based on the MIND simulation generated using GENIE and GEANT4. 
Many traversing particles, particularly hadrons, deposit energy in more than one voxel. 
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Forming clusters of adjacent voxels reduces event complexity and can improve pattern recog- 
nition in the region of the hadron shower. The clustering algorithm is invoked at the start of 
each event. The voxels of every plane in which energy has been deposited are considered in 
sequence. Where an active voxel is in contact with no other active voxel, this voxel becomes 
a cluster. If there are adjacent voxels, the voxel with the largest total deposit (at scintillator 
edge) is sought and all active voxels in the surrounding 3x3 area are considered part of 
the cluster. Adjacent deposits that do not fall into this area are considered separate. The 
cluster position is calculated independently in the x and y views as the energy-weighted sum 
of the individual voxels. One voxel, two voxel and three voxel clusters were found to have 
position resolutions of 9.4 mm, 8.0 mm and 7.2 mm, respectively. The improved resolution 
due to clusters with multiple voxels is due to the charge sharing between voxels. The clusters 
formed from the hit voxels of an event are then passed to the reconstruction algorithm. 

The separation of candidate muons from hadronic activity is achieved using two methods: 
a Kalman filter algorithm provided by RecPack [31] and a cellular automaton method (based 
on [32]), both algorithms are described in detail in [15]. The Kalman filter method requires 
a section of at least five planes where only one cluster is present in the highest z region 
of the event that is associated with particle tracks. Between 85% and 95% of (v^) CC 
interactions and ~2.5% of NC interactions fall into this category. This section is used to 
form a seed, which is projected back through the high occupancy planes using a helix model. 
Events which do not have such a section (generally high Q 2 or low neutrino energy events) are 
subject to the cellular automaton which tests a number of possible tracks to find a potential 
muon candidate. Between 5% and 13% of iv^) CC interactions and ~83% of NC are 
presented to the cellular automaton for consideration. NC events produce a candidate muon 
which is successfully fitted as such in ~60% of cases sent to the Kalman filter. Of the ~28% 

(z7 M ) CC events sent to the cellular automaton method 99% of and 45% of are 
successfully fitted. 

Compared to the method used and described in detail in [IS], an additional step has 
been added to the reconstruction method to take into account that fully-contained muons 
(particularly can have additional deposits at their endpoint due to captures on nuclei or 
due to decays. Long, well defined tracks can be rejected if there is added energy deposited 
at the muon end point, since this can be interpreted as hadronic activity and rejected by the 
Kalman filter method, thereby confusing the track finding algorithm of the cellular automa- 
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FIG. 3: Muon candidate hit purity for CC (top) and v ^ CC (bottom) interactions 
extracted using (left) Kalman filter method and (right) cellular automaton method. 

ton. Therefore, after sorting clusters into increasing z position, an additional algorithm is 
used to identify such activity and extract the track section for seeding and projection. The 
details of this algorithm can be found in [33], but it relies on identifying isolated muon-like 
hits at the end of a track and removing the high activity region in the choice of seeds to 
perform the track fit. 

The complete pattern-recognition chain using these algorithms leads to candidate purity 
(fraction of candidate hits of true muon origin) for (f m ) CC events as shown in figure [3] A 
cluster is considered to be of muon origin if greater than 80% of the raw deposits contained 
within the cluster were recorded as muon deposits. 

Fitting of the candidates proceeds using a Kalman filter to fit a helix to the candidate, 
using an initial seed estimated by a quartic fit, and then refitting any successes. Projecting 
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momentum pull p momentum (GeV/c) 

FIG. 4: Pull on the reconstructed momentum (the difference between the true and 
reconstructed momentum divided by the measured error) (left) and momentum resolution 

(right). 



successful trajectories back to the true vertex z position, the quality of the fitter can be 
estimated by comparing to the pull distribution of the reconstructed momentum, defined as 
the difference in true and reconstructed momentum, divided by the measured error in the 
momentum from the fit (see figure |4j). The error in the momentum pull is as expected, but 
the mean pull has a bias (+0.68), due to an incomplete energy-loss model in the Kalman 
filter. This small bias is taken into account in the migration matrices derived for this 
analysis (see Appendix). Further improvements to the energy- loss model within Recpack 
are being carried out and should reduce any residual bias. An empirical parametrisation of 
the momentum resolution is also shown in figure |4j which can be written as follows: 

^°' 18 + P(^- L17xir3p(Gen (1) 

Neutrino energy is generally reconstructed as the sum of the muon and hadronic energies, 
with hadronic reconstruction currently performed using a smear on the true quantities as 
described in ref. [TJ]. The reconstruction of the hadronic energy E^ad assumes a resolution 
SEhad from the MINOS CalDet testbeam E]: 

5 Eh ac i 0.55 

- 0.03. (2) 



Ehad y/Ehad 

The hadronic shower direction vector was also smeared according to the angular resolution 
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found by the Monolith test-beam [35J: 
10 - 4 i - 1 

se had = -j== © — . (3) 

V i^had ^had 

In the case of QE interactions, where there is no hadronic jet, the neutrino energy recon- 
struction was carried out using the formula: 



m N E^ + \ (m 2 N , -m^ — m 



™n — &n + \P(i \ cos v 

where $ is the angle between the muon momentum vector and the beam direction, mjy is 
the mass of the initial state nucleon, and m^i is the mass of the outgoing nucleon for the 
interactions v^ + n-^t^+p and + p — > [i + + n (see for example [SB]). The current 
algorithm only uses this formula for the case of events consisting of a single unaccompanied 
track, however, its use could be extended by selecting QE interactions using their distribution 
in i? and their event-plane occupancy among other parameters. Should the use of equation [4] 
result in a negative value for the energy, it is recalculated as the total energy of a muon with 
its reconstructed momentum. 

IV. ANALYSIS OF POTENTIAL SIGNAL AND BACKGROUND 

There are four principal sources of background to the wrong sign muon search: charge 
mis-identification of the primary muon in charged current (CC) interactions, wrong sign 
muons from hadron decay in CC events, neutral current (NC) from all species and v e CC 
events wrongly identified as CC. Typically, a charged current event has greater length 
in the beam direction than a NC or v e CC event, due to the penetrating muon. Any muons 
produced from the decay of primary interaction hadrons will tend to be less isolated from 
other hadronic activity. Additionally, the v e spectrum at a Neutrino Factory has a lower 
average energy than the spectrum which results in reduced probability of producing a 
high energy particle in the interaction. 

The previous general principles are used to define a series of offline cuts that reject the 
dominant background while maintaining good signal efficiency These can be organised in 
three categories: 1) track quality cuts; 2) charged current selection cut; and 3) kinematic 
cuts. We will describe these cuts in detail in sub-sections IV A IV B and IV C[ A summary 



of the performance of each of the cuts on signal and background will be presented in sub- 
section IIVDI and table [TTJ 



14 



g 1CT 

g 
'a 
E 

£io- 2 

o 



10" J = 



0.2 





v„ Signal 




R 




v e Background 








v NC Background 


i j 




■ v T Background 






1 : 

; j ! 


v t Contamination 




i 1 

r..i 


■ i * 'i '.; ;1 7 "'-- 


— j 



0.4 0.6 0.8 1 1.2 
Fitted Clusters/Candidate Clusters 



'io- 1 



Signal 
v e Background 

NC Background 
v, Background 
v T Contamination 




0.4 0.6 0.8 1 1.2 
Fitted Clusters/Candidate Clusters 



FIG. 5: Distribution of the proportion of clusters fitted in the trajectory for v a appearance 
(left) and appearance (right), normalised to total remaining events individually for each 

interaction type. 

A. Track quality cuts 

The quality of the reconstruction and the error on the momentum parameter of the 
Kalman filter are powerful handles in the rejection of backgrounds. We commence by im- 
posing the reconstruction criteria from the previous section to guarantee fully reconstructed 
neutrino events. We then proceed to impose a fiducial cut requiring that zl, which is the 
cluster with the lowest z in the candidate, be at least 2 m from the end of the detector 
(zl — z en d > 2000 mm), to reduce the mis-identification of candidates originating at high 
z. Additionally, a maximum value for the reconstructed muon momentum is imposed at 
16 GeV to improve energy resolution and remove backgrounds caused by very straight par- 
ticles, which confuse the fitter. 

Tracks dominated by multiple scattering or incorporating deposits made by particles not 
left by a muon can contribute significantly to backgrounds. However, these tracks will tend 
to be fitted only partially or with a larger error on the momentum variables. As such, cuts 
on these variables can be used to reduce the effect of these backgrounds. The distribution 
of the ratio of the candidate clusters which are fitted with respect to the total number of 
candidate clusters for signal and background is shown in figure [5] Accepting only those 
events in which a candidate has more than 60% of its clusters fitted reduces the background 
levels. 

Further reduction is achieved by performing a cut related to the relative error in the mo- 
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FIG. 6: Log likelihood distribution (jC q / p ) to separate wrong-sign muons for signal and 
background for v ^ (left) and (right) appearance experiments. 

mentum of the candidate muon where q is the charge of the muon and p its momentum. 
A log-likelihood distribution C q / P based on the ratio of for both signal and background 
is shown in figure |6j The signal events are selected as those with a log-likelihood parameter 
Ar/p > -0.5. 

After the preceding cuts there remain some background events which exhibit little bending 
due to the magnetic field or are reconstructed with relatively high momentum despite being 
relatively short tracks (see figure [7| right) as a result of high levels of multiple scattering. 
As can be seen in figure [7j left, removing short events in which the end point is displaced 
in the bending plane by an amount that is relatively small compared to the displacement in 
the lateral view (dispX/dispZ) effectively reduces background. Events are accepted if they 
meet the conditions described in equations [5] and [6j illustrated by the red lines in figure [7| 

d l^- > 0.18 - 0.0026 ■ iVh ; and (5) 
dispZ 

dispZ > 6000 mm or p^ < 3 • dispZ ; (6) 

where Nh is the number of clusters in the candidate, dispZ is in units of mm, and p^ in 
units of MeV/c. 

The final quality cut involves fitting to a parabola the candidate's projection onto the 
bending plane. In the current simulation a negatively-charged muon bends upwards, so that 
for a parabola defined as a + bz + cz 2 the parameter c would be positive and the charge of 
the muon is Q par = —sign(c). If the charge fitted is opposite to that found by the Kalman 
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FIG. 7: Distributions of displacement and momentum with cut levels: (top left) relative 

displacement in the bending plane to the z direction against candidate hits for signal 
events, (top right) reconstructed momentum against displacement in z for signal events 
and (bottom) as top for iv^) CC backgrounds. The red lines represent the cuts from 

equations |5j and |6j 



filter, the quality of the fit is assessed using the variable: 

if Qpar Qkal j 



QPpar 







< 


c 






< 


c 



(7) 



if Q 



par 



-Q 



kal 



where Qkal is the charge fitted by the Kalman filter fit. Defining the parameter in this way 
ensures that the cut is independent of the initial fitted charge. Events with no charge change 
(qPpar > 0.0) are accepted as signal. Additionally, those fitted badly with a charge change 
{qPpar < —1.0) are also accepted. In this way, background events which have remained in the 
sample due to local variations affecting the Kalman fitter can be removed without rejecting 
viable events in which the Kalman fitter ignored a section after a high angle scatter. The 
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FIG. 8: Distribution of the qp par variable, with the region where the parameter is < 
representing those candidates fitted with charge opposite to the initial Kalman filter for 
(left) and P M appearance experiments. The distributions are normalised to the total 
remaining events, individually, for each interaction type. 

distribution of qp par is shown in figure [8] 

B. Charged current selection 

Selection of charged currents and rejection of neutral current events is most efficiently 
performed by exploiting the propertiy that CC events tend to have greater length in z 
than NC events, since a true muon only interacts electromagnetically where a pion or kaon 
of similar momentum can interact via the strong force and will tend to stop after a shorter 
distance. Hence, the number of hits, Ihu, was used to generate Probability Density Functions 
(PDF) for charged and neutral current events (see figure [9]). One can see that the NC events 
have fewer reconstructed clusters than the equivalent CC events. For the event selection, 
candidates with greater than 150 clusters are considered signal, otherwise, the log likelihood 
rejection parameter: 

(]CC\ 

A = f || J ; (8) 



is used, which is shown in figure 10 Allowing only those candidates where the log parameter 
is C\ > 1.0 to remain in the sample ensures that the sample is pure. This analysis is similar, 
but simpler, than that employed by MINOS [37] . The effect of the CC event selection is to 
reduce the background by one order of magnitude (see table [Tl]) while having minimal effect 
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FIG. 10: Distribution of L\ likelihood ratios used to reject NC and other background 

signals. 



on the signal efficiency. 



C. Kinematic cuts 



Kinematic cuts based on the momentum and isolation of the candidate, in relation to 
the reconstructed energy of the event E rec , can be used to reduce backgrounds from hadron 
decays. The isolation of the candidate muon is described by the variable Qt = sin 2 9, 
where 9 is the angle between the muon candidate and the hadronic-jet vector. The muon 
from a true CC event is generally isolated from the hadronic jet so, on average, the Qt is 
larger for CC events than for NC events, in which a hadron associated with the hadronic 
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jet decays to a muon. Cuts based on this variable and on the reconstructed momentum 
compared to the reconstructed energy are an effective way to reduce all of the relevant beam 
related backgrounds. The distributions after the application of the preceding cuts are shown 
in figure 11, where the red lines illustrate the acceptance conditions defined in equations [9] 
andHUl 



< 5 GeV or Q t > 0.25 GeV/c and ; 



E rec < 7 GeV or 



Vp, > 0.3 -E r 



(9) 
(10) 



QE like events (see section III B ) and those events passing the conditions of equation M must 



also pass the conditions of equation 10 to remain in the data-set for the next series of cuts. 
The effect of these cuts is to reduce the background by a further order of magnitude, while 
only having a modest effect on the signal efficiency, as can be seen in table [TT| 



D. Cut summary 



In summary, after tuning the cuts described in the previous sub-sections to a test statistic, 
these were applied to independent simulated data leading to an absolute efficiency of 51% 
for selection and 62% for selection, while reducing the background to a level below 
10 -3 . A summary of all the cuts, with their effect on the signal and absolute background, 
can be found in Table |TT} The species which would be expected to contribute the greatest 
amount of background interactions for an example oscillation parameter set is also identified 
at each level. 



V. MIND RESPONSE TO THE GOLDEN CHANNEL 

Using a data-set of 3 x 10 6 events each of CC, F M CC, v e CC, V e CC and 7 x 10 6 NC 
interactions from neutrinos and anti-neutrinos generated using GENIE and tracked through 
the GEANT4 representation of MIND, the expected efficiency and background suppression 
for the reconstruction and analysis of the golden channel in MIND have been evaluated for 
both Up and z/ M appearance. Additionally, the expected level of contamination of the signal 
from other appearance oscillation channels is considered. 
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FIG. 11: Distributions of kinematic variables: (left) Reconstructed muon momentum with 
reconstructed neutrino energy for (top— >-bottom) (f m ) signal, (f m ) CC background, 
NC background, v e {v e ) CC background and (right) Q t variable (in the same order). The 
red lines represent the cuts from equations |9j and 
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TABLE II: Summary of cuts applied to select the golden channel appearance signals. The 
level of absolute efficiency and, for a 100 ktonne MIND 2000 km from the NF and 
#13 = 9.0° and 5cp — 45°, the proportion of the total non-golden channel interactions 
remaining in the sample after each cut are also shown, along with the species contributing 

the greatest number of interactions. 



Cut 


Acceptance level 


Eff. after cut 


background (xlO 3 ) 




v ii 








successful pattern rec. and fit 


0.91 


0.93 


419 (i/ e ) 


153 (Vf.NC) 


Fiducial 


zl — z en d < 2000 mm 


0.88 


0.90 


400 (y e ) 


147 (z7 M iVC) 


Max. momentum 


P M < 16 GeV 


0.85 


0.89 


158 (i/ e ) 


108 (V e ) 


Fitted proportion 


N fu /N h > 0.6 


0.81 


0.87 


74.4 (y e ) 


71.3 (i? e ) 


Track quality 


> -0-5 


0.70 


0.76 


13.6 (v e ) 


20.3 (i/ r ) 


Displacement 


dispX/dispZ > 0.18 - 0.0026^ 
dispZ > 6000 mm or < 3dispZ 


0.65 


0.72 


13.6 NC) 


10.9 (u T ) 


Quadratic fit 


qp P ar < -1.0 or qpp ar > 0.0 


0.65 


0.72 


10.3 (Uft NC) 


10.9 (i/ T ) 


CC selection 


d > 1.0 


0.63 


0.70 


2.1 NC) 


3.0 (u T ) 


Kinematic 


Erec < 5 GeV or Q t > 0.25 
E rec < 7 GeV or P M > 0.3E rec 


0.51 


0.62 


0.3 NC) 


0.9 (i/ T ) 



A. Signal efficiency and beam neutrino background suppression 

The resultant efficiencies for both polarities and the corresponding background levels 



expected for the appearance channels are summarised in Figs. [12] - 15 Numeric response 



matrices for each of the channels may be found in the Appendix. As can be seen in figure 



12 the expected level of background from CC mis-identification is around 10~ 4 , which is 
significantly below 10~ 3 at all energies for the new simulation and re-optimised analysis. 
This is also below the background levels achieved in [To], mainly due to the additional 
quality cuts. 

The background from neutral current interactions lies at or below the 10~ 3 level, with the 
high energy region exhibiting a higher level than the low-energy region due to the dominance 
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FIG. 12: Background from mis-identification of (z7 M ) CC interactions as the opposite 
polarity, (left) CC reconstructed as CC, (right) CC reconstructed as CC as a 

function of true energy. 
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FIG. 13: Background from mis-identification of NC interactions as (y^) CC 
interactions, (left) NC reconstructed as CC, (right) NC reconstructed as V ^ CC as a 

function of true energy. 

of DIS interactions. The increased particle multiplicity and greater likelihood of producing 
a penetrating pion that can mimic a primary muon are the primary reasons for this increase. 
As expected, the NC background tends to be reconstructed at low energy due to the missing 
energy (see appendix [x| . 

The background from v e (T7 e ) CC interactions is once again expected to constitute a very 
low level addition to the observed signal. This background is particularly well suppressed due 
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FIG. 14: Background from mis-identification of v e (z/ e ) CC interactions as (z/ M ) CC 
interactions, (left) v e CC reconstructed as CC, (right) v e CC reconstructed as CC as 

a function of true energy. 
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FIG. 15: Efficiency of reconstruction of (z/ M ) CC interactions, (left) CC efficiency, 
(right) z7 M CC efficiency as a function of true energy. 



to the electron shower overlapping with the hadron shower. If a particle from the hadronic 
jet decays into a wrong-sign muon, it has a lower energy and is less isolated than in the NC 
case, so the kinematic cuts suppress this background more. 

The efficiency of detection of the two polarities has a threshold lower than that seen 
in previous studies due to the presence of non-DIS interactions in the Monte Carlo sample. 



The efficiencies expected for the current analysis are shown in figure 15 



A comparison of the resultant v„ and v„ efficiency can be made with that extracted in 
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previous studies. Analyses performed in 2000 [J| and 2005 [38] assumed a 50 GeV Neutrino 
Factory, so were optimised for high energy and low values of #13. The background rejection 
achieved was at the level of 10~ 6 , but at the expense of signal efficiency, especially below 
10 GeV. A more recent analysis [9] was the first attempt at re-optimzing for a 25 GeV 
Neutrino Factory, while that in 2010 [15] was still based on a GEANT3 model, but included 
full event reconstruction and a likelihood based analysis for the first time. There exists an 
improvement in threshold for the current analysis, between 2-3 GeV, due to the inclusion 



of QE and resonance events, since these events are easier to reconstruct (see figure 20). 

The difference in efficiency between the two polarities is effectively described by the 
difference in the inelasticity of neutrino and anti-neutrino CC interactions. Neutrino DIS 
interactions with quarks have a flat distribution in the Bjorken variable: 

with Ei the scattered-lepton energy. However, anti-neutrinos interacting with quarks follow 
a oc (1 — y) 2 distribution [39|. For this reason, neutrino interactions generally involve a 
greater energy transfer to the target. The efficiencies for the two species as a function of y 



can be seen from figure 16 -(left). The shape of the efficiency curves is a consequence of the 
ratio of DIS to non-DIS events in the event samples. Neutrino and anti-neutrino efficiencies 
are very similar, showing that the Bjorken y of each event is the dominant contributor to 
the efficiency. The difference in neutrino and anti-neutrino efficiencies, when translated into 
true neutrino-energy, can be explained by the greater abundance of neutrino events at high 
y. However, since the cross section for the interaction of neutrinos is approximately twice 
that for anti-neutrinos, it is not expected that this reduced efficiency will affect the fit to 
the observed spectrum significantly. 



B. Contamination from oscillation channels containing v T or v r 

Three million events of both v T and V T interactions were generated using the GENIE 
framework [IB] and passed through the GEANT4 simulation of MIND. These events were 
then subject to the same digitization, reconstruction and analysis as the main beam back- 
grounds. Matrices were extracted describing the expected level of contamination in the 
golden channel data-set for the situation when a viable muon candidate from a v T (F r ) in- 
teraction is reconstructed as a (F^) candidate with the same and opposite charge to the 
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FIG. 16: Up CC and CC signal detection efficiency as a function of y (left) and the 
normalised distribution of all events considered in each polarity as a function of y (right). 



true primary r. As can be seen in figure 17, between 1% and 3% of the v T (y T ) interactions 
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FIG. 17: Expected level of contamination from v T (v T ) CC interactions due to the 
platinum channel, (left) v T CC reconstructed as CC, (right) V T CC reconstructed as F M 

CC as a function of true energy. 



are expected to be identified as the golden (y^) interactions. Considered properly, this 
contamination should not weaken the extraction of the oscillation parameters (see [2T]). 
Contamination from the dominant oscillation (which requires reconstruction with the oppo- 



site primary lepton charge) is expected to be below the 10 level (as shown in figure 18) 



This contamination is taken into account, but does not deteriorate the Sep fits, since the 
dominant oscillation is less sensitive to this parameter for large values of # 13 . 
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FIG. 18: Expected level of contamination from v T {p T ) CC interactions due to the 
dominant oscillation, (left) V T CC reconstructed as CC, (right) v T CC reconstructed as 

V.. CC function of true energy. 

C. Interaction expectation for 10 21 muon decays 

Using the response matrices extracted using the analysis described in the preceding sec- 
tions it is possible to make a prediction of the expected contribution to the Monte Carlo 



sample from each of the neutrino types in the beam. Figure [19] shows the expected num- 
ber of events for the best fit values of the currently measured parameters taken from [40J: 
12 = 33.5°, 23 = 45°, Am 2 21 = 7.65 x 1(T 5 eV 2 , Am 2 2 = 2.4 x 1(T 3 eV 2 for 5 CP = 45° and 
calculating matter effects using the PREM model |41| . The number of interactions were 
calculated for a 100 ktonne MIND at a distance of 2000 km from the NF for a value of 
6>i3 = 9.0°, for an integrated flux due to 10 21 decays of each polarity in the straight sections 
of the decay pipes. 



VI. STUDY OF SYSTEMATIC UNCERTAINTIES 



The efficiencies and backgrounds described above will be affected by several systematic 
effects. There will be many contributing factors including uncertainty in the determina- 
tion of the parameters used to form the cuts in the analysis, uncertainty in the exclusive 
cross-sections, uncertainty in the determination of the hadronic shower energy and direction 
resolution, and any assumptions in the representation of the detector and readout. While 
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FIG. 19: Expected interactions in a 100 ktonne MIND 2000 km from the Neutrino Factory 
for #i3 = 9.0°. Left column for appearance and right for appearance. 

exact determination of the overall systematic error in the efficiencies is complicated, an es- 
timate of the contribution of different factors can be obtained by setting certain variables 
to the extremes of their errors. 

The exclusive QE, DIS and 'other' cross-sections in the data sample could have a signifi- 
cant effect on the signal efficiencies and backgrounds. The efficiencies for the reconstruction 



of true QE and true DIS interactions are compared to the nominal efficiency in figure [20 
where the dominance of DIS interactions in the backgrounds is clear. Although experimen- 
tal data are available, confirming the presence of non-DIS interactions in the energy region 
of interest, there are significant errors in the transition regions (see for example [^2l H3]). 
These errors lead to an uncertainty in the proportion of the different types of interaction 
that can affect the efficiencies. In order to study the systematic error associated with this 
effect, events of certain types were randomly removed from the data-set and the mean effect 
quantified. As an illustration of the method, consider the contribution from QE interactions. 
Taking the binned errors on the cross-section measurements from |l2j H3] , a run to reduce 
the QE content would exclude a proportion of events in a bin so that instead of contribut- 
ing the proportion where Nqe and N tot are the number of QE interactions and the 
total number of interactions in the bin of interest, it would instead contribute: UQE ^ QE 
where oqe is the proportional error on the QE cross section for the bin. Since the data-set is 
finite and an actual increase in the number of QE interactions is not possible, the equivalent 
run to increase the QE contribution reduces the contribution of the "rest" by an amount 
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FIG. 20: Efficiencies for a pure DIS sample compared to the nominal case, (top) Signal 
efficiency, (second line) (f m ) CC background, (third line) NC background and (bottom) 
v e {y e ) CC background, appearance on the left and appearance on the right. 
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calculated to give the corresponding proportional increase in QE interactions: 
Nqe + o'qeNqe _ N QE 



(12) 



N t ot + vqeNqe Ntot - eN rest ' 
where N rest is the total number of non-QE interactions in the bin and e is the required 
proportional reduction in the 'rest' to simulate an appropriate increase in QE. Solving for e 
yields the required reduction: 

e = ^^. (13) 

The la systematic error can be estimated as the mean difference between the nominal 
efficiency and the increase due to a higher QE proportion or decrease due to exclusion. 
The errors in the true and F M efficiencies extracted using this method that varies the 



contribution of QE, \tx and other non-DIS interactions are shown in figure 21 Errors for 
l7r resonant reactions are estimated to be ~20% below 5 GeV (as measured by the K2K 
near detector jH]) and at 30% above. Due to the large uncertainty, both theoretically and 
experimentally, on the models describing other resonances, coherent, diffractive and elastic 
processes, a very conservative error of 50% is taken when varying the contribution of the 



"others". As can be seen in figure 21, the systematic effect is less than 1% for neutrino 
energies between 3 GeV and 10 GeV increasing to 4% near 1 GeV, with increased QE 
and l7r interactions generally increasing the efficiency and increased contribution of the 
"other" interactions having the effect of decreasing efficiency. This last result is likely to be 
predominantly due to resonances producing multiple tracks. The effect on backgrounds is 



expected to be minimal, as was also shown in figure 20 At the time of a Neutrino Factory, 
the cross-section uncertainties should be much smaller than the ones assumed here, so we 
expect the systematic error to be below 1% for all energies. 

Another important source of systematic uncertainty is due to the error in the muon mo- 
mentum and the hadron shower energy used to reconstruct the total neutrino energy. The 
muon is fully reconstructed, but the hadron shower reconstruction is performed assuming 
a parametrisation to the hadron energy and angular resolution. The particular choice of 
this parametrisation introduces a systematic error on the resulting signal and background 
efficiencies that should be studied. Taking a 6% error as quoted for the energy scale uncer- 
tainty assumed by the MINOS collaboration [7] and varying the constants of the energy and 



direction smears by this amount, it can be seen (blue bands in figure 22) that, to this level, 
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FIG. 21: Calculated error on signal efficiencies on increasing (top) and decreasing 
(bottom) the proportion of non-DIS interactions in the data-set. (left) Errors on true 
energy CC efficiency and (right) errors on true energy V ^ CC efficiency 



the hadronic resolutions have little effect on the true neutrino-energy efficiencies. However, 
the hadronic direction resolution is likely to have far greater uncertainty and would be very 



sensitive to noise in the readout electronics. Also shown in figure [22] are the efficiencies 
when the hadronic energy resolution parameters are 6% larger but with a 50% increase in 
the angular resolution parameters. We expect in a real detector to measure the hadronic 
angular resolution with a precision of better than 50%, even though we cannot quantify this 
precision yet. However, even at this level, the observed difference in efficiency is only at 
the level of 1% above 7-8 GeV. A combination of the exclusive cross-sections and hadronic 
energy uncertainties implies a total systematic uncertainty for the measurement efficiency 
of order 1% over the neutrino energy range above 2 GeV. 
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FIG. 22: Variation of signal efficiency (top) and NC backgrounds (bottom) due to a 6% 
variation in the hadron shower energy and direction resolution and a more pessimistic 50% 
reduction in angular resolution (focused on region of greatest variation). 

VII. COMPARISON BETWEEN EVENT GENERATORS 



The above analysis uses the GENIE event generator to simulate the interaction of events 
with matter in the detector. This generator was assumed to bring this effort in line with 
current experiments, such as MINOS, where good agreement with data has been achieved. It 
is useful to compare to a previous version of this analysis with the NUANCE event generator 



20] . A comparison of the neutrino charge current detection efficiencies appears in figure 23 



This shows that the GENIE derived analysis produces smaller positive identification of 
charge current events. This loss of performance is linear with respect to neutrino energy for 
energies greater than 5 GeV, so that the analysis of the GENIE simulation is 20% less efficient 
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(b) GENIE and NUANCE P^CC efficiency 
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FIG. 23: The Golden channel analysis efficiencies from simulations of charge current 
neutrino events in a MIND detector generated using GENIE and NUANCE neutrino 
generators. The differences between the efficiencies as fractions of the GENIE efficiency are 

also shown. 



than that of the NUANCE simulation. This difference is partially ascribed to differences in 
the parton distribution functions used by each generator and should not be interpreted as a 
systematic error of the analysis, since the measured event rates in a future experiment will 
be cross-checked against the simulations, so the efficiencies will be determined much more 
accurately than the difference between generators. We assume that GENIE, which has been 
benchmarked against recent neutrino experiments, serves as a more realistic estimator for 
the efficiency of the analysis. 
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VIII. MIND SENSITIVITY 



The Neutrino Factory is required to measure the CP-violating phase 5cp simulaneously 
with 9i3 while removing ambiguity caused by degenerate solutions. Extracting the oscillation 
parameters from the observed signal at the far detectors requires the accurate prediction 
of the expected flux without oscillation, which is used along with the calculated oscillation 
probabilities to fit the observed signal for the best value of #13 and 5c p- Due to the large 
distance to the far detectors, the flux spectra expected from the decay rings of the Neutrino 
Factory are accurately approximated by the flux from a point source of muons travelling with 
appropriate Lorentz boost in the direction of the detector. Using this flux or, alternatively, 
a projection of the spectra observed in the near detector, the number of true interactions 
expected in MIND as a function of energy can be calculated for each value of #13 and S. 
These spectra can then be multiplied by the response matrices shown in the Appendix to 
calculate the observed golden channel interaction spectrum expected for some hypothetical 
values of the oscillation parameters. 

A Neutrino Factory storing muons of energy 10 GeV is assumed, of which 5.0xl0 20 per 
year of each species decay in the straight sections pointing towards the MIND far detector 
of 100 ktonnes mass placed at a distance of 2000 km from the facility. 



A. The NuTS framework 

The Neutrino tool suite (NuTS) was developed for the studies presented in [iBT - HT] . It 
provides a framework for the generation of appropriate fluxes for different neutrino accelera- 
tor facilities along with the necessary infrastructure to calculate the true neutrino oscillation 
probabilities for all channels. In addition, using the parametrisation of the total interaction 



spectra calculated in section HI A the expected number of events in a given energy bin 
can be calculated. Using this framework and the response matrices extracted for MIND, 
simulated data for an experiment can be generated as 



Data± = smear M\ ig N% + £ M*N$ , (14) 



\ k / 

for each polarity and detector baseline of interest, where M\ ig is the response matrix for 
MIND for a particular signal channel % (stored fi + or //~), N l J g is the 100% efficiency inter- 
action spectrum in true v energy bins for a channel i at a detector baseline j (in this case, 
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there is only one baseline at 2000 km), M^ g is the response matrix for a background k (mis- 
identification of CC interactions from other neutrino species or from NC) to the appearance 
channel i, N^ g is the 100% expectation spectrum for a background k to an appearance 
signal i at a detector baseline j and these expected values are used to calculate an observed 
number of interactions following a Poisson distribution. 



B. Fitting for #13 and 5c p simultaneously 

Due to the correlation between #13 and 5c p, a simultaneous fit is necessary. Defining a 
grid of 6*i3 and 5c p values the \ 2 °f a fit to Dataf im can be calculated using the function 



X 2 = E { 2 x E (■{,■•:, y.JOv.', S C p) - n e +J + n%j log ( : 



AjXjN* tj (8 13 ,8 C p) 



+A j Ni J (9 vl , 5 C p) - nt d + ntj log y AjN1Jei3: s Cp) 
, (Aj - I) 2 , (*j ~ I) 2 



0"A &x 



(15) 



where nf • is the simulated data (Datdf im ) for an energy bin e, Nfj{9\^5cp) is the pre- 
dicted spectrum for the values of 6*13 and 5c p represented by the grid point (calculated as 



in equation 14 but without a smear) and j represents the baseline as in equation 14 The 
uncertainty in the expected number of interactions and expected ratio in interactions be- 
tween neutrinos and antineutrinos are represented by the additional free parameters Aj and 
Xj respectively and their corresponding errors. 

We made two assumptions regarding the overall event normalisation: one assumes a 
conservative error of a a = 0.025 and the other assumes a more optimistic (but realistic) 
assumption for a neutrino factory of oa — 0.01. The uncertainty in the ratio of cross 
sections between neutrinos and antineutrinos is maintained fixed at a x = 0.01, which is 
the level to which a near detector would seek to measure the interaction cross-sections at 
the time of a neutrino factory. The minimisation of the parameters A and x is performed 
analytically for each predicted dataset to leading order. The contours at Xmin + 9 represent 
approximately the 3<r level of understanding and those at Xmin+25 represent 5a. In such a fit, 
the experimentally determined oscillation parameters jlQ] are considered fixed. While there 
would be some systematic error associated with uncertainty in these parameters, systematics 
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FIG. 24: Examples of la, 3a and 5a fits to simulated data for the 10 GeV Neutrino 

Factory, for sin 2 26 13 = 0.096. 



from the normalisation and cross-section uncertainties are expected to dominate. Some 
examples for such fits, assuming sin 2 2#i3 = 0.096, are shown in figure 24 



C. Sensitivity to the mass hierarchy 

Distinguishing between the two possible mass hierarchies is important for the understand- 
ing of the neutrino mass sector. The sensitivity to the true mass hierarchy is defined here as 
when the true sign of Am 2 3 can be distinguished from the opposite sign to the appropriate 
na level, that is 

xLn(- Am L) - xLn( Am is) > n2 - ( 16 ) 

Sensitivity to the mass hierarchy up to 5a covers the whole region down to sin 2 2#i3 < 10~ 4 , 
so the Neutrino Factory would be able to cover easily the current value of sin 2 2#i 3 ~ 0.1. 



D. Sensitivity to Sep 

Sensitivity to the measurement of 5c p is defined as when the difference between the 
minimum y 2 with respect to the minimum obtained by fitting with CP conserving cases is 
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FIG. 25: Sep 3a, 5a and 10a measurements (top) and 5cp 5a coverage to measure CP 
violation (bottom) as a function of sin 2 29\3 for true normal hierarchy (left) and true 
inverted hierarchy (right). The vertical lines represent the range of possible values of 

sin 2 20i3. 

greater than the appropriate na level, that is 

rmn(x 2 (6 CP = 0), X \^cp = 180), X 2 (S CP = -180)) - xLn > n 2 . (17) 
The 5c p measurement and 5c p coverage plots to measure CP violation by the 10 GeV 



Neutrino Factory are shown in figure 25 for the normal mass hierarchy (left) and for the 
inverted mass hierarchy (right). The CP coverage to measure 5cp at the ha level for both 
normal and inverted mass hierarchy is ~85% in the range of the currently measured values 
of sin 2 2#i3. These results improve those presented by the Interim Design Report [20] . 

The accuracy achieved in the measurement of 5c p is an increasingly important parameter 
to determine the performance of a facility [H]. The la error A5cp for #13 = 9° is shown 



in figure 26 (left), under two assumptions of the overall normalisation and cross-section 
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FIG. 26: Left: Precision in the value of 5cp as a function of Sep under two different 
assumptions: top line (3.0% normalisation and 2.5% cross-section error), bottom line (1% 
normalisation and cross-section error). Right: 5c p fraction coverage that can be achieved 
above each value of A5cp for 3.0% and 2.5% errors (right curve) and for 1% and 1% errors 

(left curve). 

systematic errors: (<JA,cr x ) = (3.0%, 2. 5%) and (1.0%, 1.0%). Depending on the value of 
5qp and the level of systematic error, the accuracy in A5cp is between 2.5% and 5%. 



Figure 26 (right) shows the 6 fraction coverage that can be achieved above the value of AScp 



determined at a Neutrino Factory for each of the two assumptions about the systematic error. 



IX. CONCLUSIONS 

A new GEANT4 simulation of the Magnetised Iron Neutrino Detector (MIND), using 
GENIE as the neutrino event generator, has been developed to determine the performance 
of MIND at a Neutrino Factory. Considering the spectrum of neutrino interactions in the 
energy region 0-25 GeV produced by the GENIE generator, this simulation has been used 
to study the efficiency and the background rejection of MIND. A detector of 100 ktonne 
at 2000 km was used to determine the expected sensitivity to 5c p and the mass hierarchy 
at a Neutrino Factory for a value of sin 2 9 13 ~ 0.1, as determined recently by reactor and 
accelerator neutrino oscillation experiments. 

The proportions of quasi-elastic, single pion production and deep inelastic events obtained 
from GENIE have been benchmarked with experiments and a parametrisation of the total 
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interaction spectrum, in agreement with data, was used for the simulation. Digitisation of 
the events tracked through the simulation assumed read-out of the scintillator using WLS 
fibres and electronics with 30% QE and a standard deviation on signal response of 6%. 
Both assumptions should be achievable and could, perhaps, turn out to be conservative if 
the current trend in photon detectors and electronics performance continues. These events 
have been used to study the efficiency and background suppression of MIND in a search 
for wrong sign muons at a Neutrino Factory storing both muon polarities. A re-optimised 
reconstruction algorithm and a new analysis applied to data-sets comprising several million 
v fi ipn) an d v e (u e ) CC and NC events, resulted in response matrices describing the expected 
response of the detector. MIND showed an efficiency plateau from 5 GeV in true neutrino 
energy at ~60% for and ~70% for and thresholds at ~2 GeV. All beam inherent 
backgrounds were simultaneously suppressed at a level around 1CT 4 . The difference in 
efficiencies for the two polarities has been studied and found to be predominantly due to the 
difference in the inelasticity spectrum expected for neutrino and anti-neutrino interactions. 
We expect the dominant systematic errors to be due to the hadronic energy resolution and 
due to the different neutrino interaction types. Each contribute an expected systematic error 
on the signal efficiency of around 1%. 

The response of MIND and a parametrisation of the total neutrino and antineutrino cross- 
sections, calculated from the output of the GENIE event generator, were used to study the 
sensitivity of two MIND detectors at the Neutrino Factory. The obtained measurement of 
key parameters of the PMNS matrix indicate that such an experiment would determine the 
mass hierarchy, irrespective of the value of Sep, at a level better than 5a, could perform 
measurements of Sep with an accuracy between 3% and 5% and would have an 85% Sep 
coverage for the currently preferred value of 9 13 ~ 9.0°. 
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X. APPENDIX 

This appendix summarizes the response matrices of signal (wrong sign and P M appear- 
ance) and all backgrounds, [p^ and CC, u e and u e CC, z/ r and u T CC, and NC) in bins 
of true and reconstructed neutrino energy relevant to an oscillation analysis. Each entry in 
the table is the survival probability for each species. In all tables, columns represent the 
true neutrino energy in GeV and rows the reconstructed energy, also in GeV. The overflow 
bin in reconstructed energy represents all events with a reconstructed energy greater than 
the known maximum. 

A. appearance 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 





1.392 


3.940 


2.052 


2.402 


1.323 


0.3669 


0.2460 


0.3710 


0.1218 


1.0-2.0 





222.0 


243.6 


45.87 


28.82 


21.65 


9.660 


6.519 


1.608 


2.679 


2.0-3.0 





149.4 


1238 


568.5 


144.5 


68.80 


34.73 


22.63 


14.10 


7.551 


3.0-4.0 





14.13 


807.9 


1790 


813.2 


270.8 


110.4 


60.52 


41.31 


26.79 


4.0-5.0 





1.392 


162.0 


1340 


1848 


859.0 


318.7 


146.9 


81.87 


57.24 


5.0-6.0 





0.1070 


14.64 


375.5 


1516 


1695 


805.7 


339.7 


152.2 


83.54 


6.0-7.0 





0.1070 


1.914 


62.40 


593.3 


1630 


1616 


841.6 


369.3 


178.8 


7.0-8.0 








0.3378 


7.725 


123.1 


744.8 


1563 


1462 


783.8 


347.1 


8.0-9.0 











0.7242 


18.37 


200.4 


836.3 


1547 


1360 


771.6 


9.0-10.0 








0.1126 


0.3621 


2.402 


35.24 


277.6 


891.3 


1454 


1289 


10.0-11.0 














0.6004 


8.059 


83.52 


474.8 


1520 


3014 



TABLE III: Golden channel appearance signal efficiency All values xlO 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 
































2.0-3.0 

















0.1203 


0.1218 








0.1198 


3.0-4.0 

















0.2406 








0.1220 


0.1198 


4.0-5.0 




















0.2437 


0.2439 





0.1198 


5.0-6.0 


























0.1220 





6.0-7.0 























0.1220 


0.1220 


0.2396 


7.0-8.0 

















0.1203 





0.1220 





0.2396 


8.0-9.0 





























0.1198 


9.0-10.0 




















0.1218 











10.0-11.0 





























0.1198 



TABLE IV: \x background from charge mis-identified CC events All values x 10 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 
































2.0-3.0 
































3.0-4.0 














0.0596 

















4.0-5.0 











0.0596 








0.0614 











5.0-6.0 





























0.0610 


6.0-7.0 
































7.0-8.0 
































8.0-9.0 
































9.0-10.0 
































10.0-11.0 





























0.0610 



TABLE V: /i background from u e CC events. All values xlO 4 . 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 











0.1729 


0.0581 


0.3538 


0.0602 


0.1221 





0.0612 


2.0-3.0 








0.1068 





0.3488 


0.2359 


0.1807 


0.4272 


0.3680 


0.1224 


3.0-4.0 











0.2305 


0.8138 


0.8844 


0.8431 


0.6713 


0.6133 


0.8566 


4.0-5.0 











0.1152 


0.5232 


0.7076 


1.144 


1.098 


0.9812 


0.9178 


5.0-6.0 











0.0576 


0.0581 


0.5307 


0.3011 


0.5492 


0.6746 


0.6119 


6.0-7.0 

















0.1769 


0.3011 


0.4272 


0.5519 


0.4895 


7.0-8.0 




















0.0602 


0.4272 


0.2453 


0.4283 


8.0-9.0 























0.0610 


0.1840 


0.4283 


9.0-10.0 





























0.1836 


10.0-11.0 





























0.1836 



TABLE VI: \T background from NC events. All values xl0~ 4 . 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 














0.5833 








0.3747 





0.1426 


1.0-2.0 











14.70 


8.750 


8.569 


12.33 


8.992 


9.264 


6.845 


2.0-3.0 











49.97 


51.33 


54.27 


65.98 


67.44 


52.33 


45.35 


3.0-4.0 











38.21 


51.92 


58.08 


71.91 


76.06 


78.66 


79.71 


4.0-5.0 











11.76 


18.08 


38.08 


47.49 


57.32 


61.27 


70.16 


5.0-6.0 














5.833 


16.19 


28.77 


32.41 


33.32 


42.78 


6.0-7.0 

















3.808 


12.33 


17.23 


18.69 


25.67 


7.0-8.0 

















1.587 


4.109 


10.30 


10.89 


13.69 


8.0-9.0 




















1.142 


2.810 


5.526 


7.986 


9.0-10.0 




















0.2283 


0.3747 


2.113 


4.278 


10.0-11.0 


























1.788 


2.139 



TABLE VII: n reconstructed from v T CC events. All values xlO 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 
































2.0-3.0 
































3.0-4.0 























0.1637 


0.4287 





4.0-5.0 























0.1637 


0.2858 





5.0-6.0 
































6.0-7.0 
































7.0-8.0 
































8.0-9.0 
































9.0-10.0 
































10.0-11.0 

































TABLE VIII: \T background from v T CC events. All values xlO" 4 . 



B. appearance 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 





4.606 


22.28 


14.26 


9.296 


3.729 


1.340 


0.3659 





0.7188 


1.0-2.0 





85.70 


262.9 


93.41 


54.43 


26.11 


18.15 


9.025 


4.757 


1.318 


2.0-3.0 





84.87 


1565 


656.4 


116.1 


51.85 


39.48 


28.05 


20.61 


12.94 


3.0-4.0 





5.304 


990.1 


2364 


795.4 


164.5 


63.85 


49.15 


46.11 


38.57 


4.0-5.0 





0.1396 


161.9 


1904 


2513 


880.2 


226.6 


89.76 


66.35 


54.63 


5.0-6.0 








10.52 


477.7 


2145 


2218 


851.9 


246.6 


97.94 


64.45 


6.0-7.0 











60.02 


801.6 


2221 


2032 


842.4 


283.0 


111.5 


7.0-8.0 








0.1238 


3.751 


152.5 


1022 


2181 


1896 


840.2 


307.4 


8.0-9.0 











0.6252 


19.33 


263.2 


1169 


2108 


1738 


832.3 


9.0-10.0 











0.1250 


1.223 


47.28 


355.9 


1230 


1996 


1633 


10.0-11.0 











0.1250 


0.3670 


8.301 


98.33 


603.4 


2024 


4033 



TABLE IX: Golden channel appearance signal efficiency. All values xlO 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 
































2.0-3.0 
































3.0-4.0 








0.3378 


0.6035 





0.2406 


0.1223 


0.1230 








4.0-5.0 








0.2252 


0.8449 


0.3602 


0.2406 











0.2436 


5.0-6.0 











0.1207 


0.2402 





0.1223 


0.3690 


0.1237 


0.1218 


6.0-7.0 

















0.1203 


0.2446 





0.1237 





7.0-8.0 

















0.2406 


0.1223 


0.1230 








8.0-9.0 











0.1207 




















9.0-10.0 























0.1230 


0.1237 


0.1218 


10.0-11.0 

















0.1203 


0.2446 


0.2460 





0.3654 



TABLE X: fi + background from charge mis-identified CC events All values x 10 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 
































2.0-3.0 

















0.0605 


0.0608 











3.0-4.0 








0.0613 





0.0609 








0.1221 


0.1215 


0.0602 


4.0-5.0 











0.0627 


0.1218 


0.1211 


0.2431 


0.1832 


0.1215 


0.1204 


5.0-6.0 























0.0611 








6.0-7.0 














0.0609 








0.1221 








7.0-8.0 
































8.0-9.0 





























0.0602 


9.0-10.0 























0.0611 


0.0608 


0.0602 


10.0-11.0 

































TABLE XI: fi + background from v e CC events. All values xlO 4 . 
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0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 

















0.0586 





0.0607 








2.0-3.0 











0.1167 





0.1173 





0.0607 





0.0606 


3.0-4.0 

















0.1759 


0.0597 


0.1213 


0.1210 


0.2422 


4.0-5.0 

















0.0586 


0.1194 


0.2427 


0.1814 


0.0606 


5.0-6.0 




















0.0597 


0.0607 


0.1210 


0.1817 


6.0-7.0 























0.0607 


0.1814 


0.0606 


7.0-8.0 


























0.0605 





8.0-9.0 
































9.0-10.0 





























0.0606 


10.0-11.0 

































TABLE XII: fj+ background from NC events. All values xl(T 4 . 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 











15.53 


1.498 


2.763 


3.066 


2.128 


1.857 


2.020 


1.0-2.0 














22.47 


33.16 


31.88 


27.66 


25.29 


23.99 


2.0-3.0 











62.11 


64.41 


77.99 


74.19 


75.28 


67.44 


59.34 


3.0-4.0 











54.35 


63.67 


68.47 


66.62 


61.70 


73.59 


69.82 


4.0-5.0 











15.53 


28.46 


42.99 


54.36 


55.15 


55.15 


53.91 


5.0-6.0 














5.243 


19.65 


36.17 


36.99 


41.29 


43.43 


6.0-7.0 














3.745 


5.220 


11.44 


20.62 


27.01 


26.39 


7.0-8.0 

















0.6141 


4.292 


12.27 


16.72 


18.94 


8.0-9.0 




















1.635 


2.619 


5.430 


10.48 


9.0-10.0 




















0.4087 


1.964 


3.572 


5.555 


10.0-11.0 

















0.3070 


0.2044 


0.4910 


1.286 


5.176 



TABLE XIII: fi + reconstructed from v T CC events. All values x 10 





0.0-1.0 


1.0-2.0 


2.0-3.0 


3.0-4.0 


4.0-5.0 


5.0-6.0 


6.0-7.0 


7.0-8.0 


8.0-9.0 


9.0-10.0 


0.0-1.0 
































1.0-2.0 




















0.6849 


0.1873 


0.6501 


0.2852 


2.0-3.0 

















0.3174 


0.2283 


0.9367 


0.4876 


0.9982 


3.0-4.0 














0.5833 





0.2283 


0.3747 


0.9751 


1.426 


4.0-5.0 























0.3747 


0.4876 


0.7130 


5.0-6.0 


























0.3250 


0.8556 


6.0-7.0 





























0.1426 


7.0-8.0 
































8.0-9.0 
































9.0-10.0 
































10.0-11.0 

































TABLE XIV: /i+ background from v T CC events. All values xlO" 4 . 



